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We present first measurements of the pseudorapidity and azimuth (ij, <f>) bin-size dependence 
of event-wise mean transverse momentum (p t ) fluctuations for Au-Au collisions at ^/snn = 200 
GeV. We invert that dependence to obtain p t autocorrelations on differences (»ja,0a) interpreted 
to represent velocity/temperature distributions on (?7, </>). The general form of the autocorrelations 
suggests that the basic correlation mechanism is parton fragmentation. The autocorrelations 
vary rapidly with collision centrality, which suggests that fragmentation is strongly modified by a 
dissipative medium in the more central Au-Au collisions relative to peripheral or p-p collisions. 
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Central Au-Au collisions at RHIC may generate 
a color-deconfined medium (quark-gluon plasma or 
QGP) PJ. Some theoretical descriptions predict abun- 
dant low-pt gluon production in the early stages of high- 
energy nuclear collisions, with rapid parton thcrmaliza- 
tion as the source of that medium [2j,|3J,|4|. Particle yields, 
spectra and high-p t correlations from Au-Au collisions at 
^snn = 130 and 200 GeV provide tantalizing evidence 
that a colored medium is produced [a, IS Li B- Nonstatis- 
tical fluctuations of event- wise mean-p t (p t ) ^[UE3 should 
help to determine the properties of that medium. A re- 
cent measurement of excess (pt) fluctuations in Au-Au 
collisions at 130 GeV revealed a large excess of fluctu- 
ations compared to independent-particle pt production. 
The measurement was obtained at a single scale (bin 
size) — the STAR detector acceptance on (77, (f>) for that 
analysis. Excess (p t ) fluctuations studied with Monte 
Carlo simulations have been attributed to low-p t parton 
fragments (minijets) Measurements of (p t ) fluctu- 

ations could help to illuminate the role of minijets in 
nuclear collisions. 

In this paper we report the first measurements of the 
scale dependence of (pt) fluctuations. Moreover, by in- 
version of the scale-dependent {pt) variance distribution 
we obtain p t autocorrelations, projections of two-particle 
distributions on momentum difference variables (r]A, 4>a), 
where e.g., t]a = Vi ~~ V2 LL2i- We compare the result- 
ing p t correlation patterns to known azimuthal corre- 
lations (e.g., elliptic flow) and jet angular correlations. 
We consider the possibility that minijets, as local veloc- 
ity correlations, provide a dominating contribution to p t 
correlations and quantify centrality dependencies which 
may describe in-medium modification of jet correlations. 
This analysis is based on ^/snn = 200 GeV Au-Au colli- 
sions observed with the STAR detector at the Relativistic 
Heavy Ion Collider (RHIC). 

In each heavy ion collision, and within some region 
on (77, 4>) called a bin, a number of individual particle 
p t s is sampled from local pt spectra. Local spectrum 
properties may deviate from the event-averaged pt spec- 
trum differently at each point on (77, cf) and differently in 
each event 13] . The bin-size (scale) dependence of excess 
event- wise (p t ) fluctuations measured by variance differ- 
ence A<7p t . n (5?7, S(f>) reflects the correlation structure of 
the local pt spectrum properties . Certain aspects of 
the correlation structure can be accessed when that scale 
dependence is inverted to obtain pt autocorrelations ^3] : 
those aspects which depend on relative separation of pairs 
of points but not on absolute position on (r],4>). The p t 
autocorrelations for Au-Au collisions over a range of cen- 
tralities, their structure and interpretation, are the main 
subjects of this paper. The next three paragraphs define 



the (pt) fluctuation measure and outline the derivation of 
the integral equation which connects its scale variation 
to the corresponding autocorrelation distribution. Those 
paragraphs may be omitted in a first reading. 

In this analysis the detector acceptance is divided into 
macrobins with scales (5rj,8<j)). Each macrobin (scales 
represented by Sx for brevity) contains in each event some 
integrated particle multiplicity n(Sx) and total pt(Sx) 
(scalar sum over particles in the bin). Rather than fluc- 
tuations of ratio (pt) = pt/n (a source of systematic er- 
ror) we study fluctuations of difference (p t — npt)j\ffi. 
The scale-dependent per-particle pt variance is defined by 
o~p t . n {Sx) = (p t (Sx) — n(Sx) pt) 2 /n(8x), where pt is the 
inclusive mean particle p t , n is the mean bin multiplicity, 
Pt '■ n reads pt given multiplicity n, and the overline rep- 
resents an average over all macrobins in all events [lOj. 
The small-scale limit (n = 1) of (Jp t . n {8x) is er| t , the in- 
clusive single-particle variance. The variance difference is 
then defined as Ao-p t . n (5x) = o-p t . n (Sx) — cr| t . Variation 
of AcTp t . n on scales (5r],S(j)) corresponds to an integral 
equation which can be inverted to obtain p t autocorrela- 
tions on difference variables (?7a,0a), which compactly 
represent two-particle p t correlations on (77, <p) [15j and 
permit direct interpretation of (p t ) fluctuations in terms 
of physical mechanisms. 

The autocorrelation distribution is a powerful tool for 
accessing two-particle correlations under certain condi- 
tions well satisfied in relativistic nuclear collisions jl5| . 
An autocorrelation compares a distribution f(x) to it- 
self. It is effectively a projection by averaging of prod- 
uct distribution f(xi) ■ /(a^) on (xi,x 2 ) onto the differ- 
ence variable xa = x\ — x 2 . In this analysis we obtain 
the autocorrelation of the pt distribution on 2D space 
(77, <t>). Autocorrelations can be determined by pair count- 
ing [lHl ITfil ] , or by inverting fluctuation scale dependence 
to form density ratios following the procedure in [l2T | first 
implemented in |17|. Here we use the latter method. 

We can relate variance measurements to autocorre- 
lations in the following way. If a space x is par- 
titioned into microhms of fixed size e x , combined to 
form macrobins of variable size Sx, the macrobin con- 
tents in a 2 t . n (5x) can be expressed as microbin sums, 
e.g., p t {5x) = J2 a PtA e x)i a being a microbin index. 
We can then express variance o~ 2 t . n (5x) as a double 
sum over microbin indices (a, 6) on (x\,X2) of terms 
(pt — npt) a {pt — np t )b, which measure the covariance be- 
tween bins a and b on x of pt fluctuations relative to 
npt . As shown in [T2I Il7| we can rearrange the dou- 
ble sum into an outer sum over index k on difference 
variable xa = x\ — X2 {e.g., r\\ — 772, with microbin in- 
dex k) and an inner sum over microbins on sum variable 
xi + X2- The inner sum is pt difference autocorrelation 
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AAk(pt ■ n) ('difference' referring to pt — np t ). If self 
pairs are excluded from the microbin sums the pt differ- 
ence autocorrelation corresponds to variance difference 
Aa pt . n (Sx). We define reference number autocorrelation 
Ak,ref (n) as the mean pair number n 2 , in the k th microbin 
on obtained by averaging products of mean particle 
numbers n a fib along the k th diagonal of (xx, X2), that is, 
with a — b = k. That reference is approximately equiv- 
alent to the mixed-pair reference autocorrelation which 
would be obtained by direct pair counting Ak tre f{n) 
is not obtained explicitly in this analysis, is instead an 
implicit part of the density ratio defined below and in- 
ferred by fluctuation inversion. Autocorrelation densities 
p(xa), defined e.g. by AA k (p t : n) = e 2 x Ap(p t : n;ke x ) 
and Ak ir ef(n) = e 2 p re f(n; ke x ), are independent of mi- 
crobin size. The required per-particle autocorrelation 
measure corresponding to Ao^ t . n (fa) is density ratio 

Ap(p t : n)/yjpref{n) = AA{p t : n)/e x yj 1 A ref {n) [units 
(GeV/c) 2 ], which estimates the p t covariance per parti- 
cle for a given separation on (n,(j)), averaged over the 
acceptance 0]. Within an 0(1) constant factor such 
density ratios have the form of Pearson's correlation co- 
efficient |19| : the average covariance for all pairs of bins 
with a given separation on (77, <f>) relative to the geometric 
mean of Poisson number variances for those bin pairs. 

For this 2D scaling analysis we generalize Sx — > 
(Srj, 5(f>) to obtain the per-particle conditional p t variance 
difference (also defining difference factor Aa Pt:n 0[) as 
the 2D discrete integral equation 

^o-l t .n{mse rn n s e (t> ) = 2a Pt Aa Pt . n (m s e v , n s e^) (1) 



. K Ap{p t 



n-^ke^Jccf 



k,l=l 



with kernel K, 



\/pref(n;ke n ,le^) 

(m s - k+l/2)/ms ■ {ns - 
I + 1/2) /ns- That integral equation can be inverted to 
obtain autocorrelation density ratio Ap/^/p re f as a per- 
particle pt correlation measure on (t)a, 4>a) from the scale 
dependence of (pt) fluctuations represented by variance 
difference Aa 2 pt , n {5n,8cj)) 0H3. 

Data for this analysis were obtained with the STAR 
detector [2lJ using a 0.5 T uniform magnetic field par- 
allel to the beam axis. Event triggering and charged- 
particle measurements with the time projection chamber 
(TPC) are described in [2l]]. Track definitions, tracking 
efficiency and background corrections, event and track 
quality cuts and primary-particle definition are similar to 
those described in [l(J, |22| for 130 GeV. While there are 
quantitative differences between the two energies and de- 
tector configurations, the better quality of full-magnetic- 
field tracking at 200 GeV tends to more than offset the 
effect of larger track densities there compared to the half- 
field tracking at 130 GeV. The difference between back- 
grounds is a few percent of the total track yield (larger 
for 200 GeV) and is included in the corrections. Tracks 
were accepted with pseudorapidity in the range < 1, 
transverse momentum in the range pt € [0.15, 2] GeV/c 



and 27r azimuth, defining the detector acceptance for this 
analysis. Particles identification was not implemented. 
Eleven centrality classes were defined as fractions of a tot 
(nine equal fractions from 90% to 10%, the top 10% being 
further divided in half). The centralities specified below, 
rounded to the nearest 5%, are within 2% of the defined 
values. Centralities were determined using the uncor- 
rected number N of charged particles in \rf\ < 1 [23|. 




FIG. 1: Left panels: Aa Pt:n (GeV/c) 2 distributions on scale 
(577, 8(f)) for three centrality bins: 80-90% of total cross section 
(top panel), 45-55% of total cross section (middle panel) and 
0-5% of total cross section (bottom panel). Right panels: Cor- 
responding autocorrelations on difference variables (?7a,0a)- 

Fig- (left panels) shows the scale dependence of vari- 
ance difference Aa^ t . n (8T), £</5) in Eq. ^ f° r three of the 
eleven centralities analyzed. The scale axes are divided 
into microbins: 16 on pseudorapidity scale Sn and 24 on 
azimuth scale 6cj>. Variance differences typically increase 
monotonically with Sn but have more complex behav- 
ior on 5(f>. Measurements of difference factor A<r Pt:n at 
y/SNN = 130 GeV reported in 0] correspond to the sin- 
gle point at the STAR acceptance scale (2, 2ir) for each 
centrality. To access the underlying dynamics we extract 
the corresponding autocorrelation distributions. Fig. ^ 
(right panels) shows 2D autocorrelations (by construc- 
tion symmetric about »7a,<M = 0) inferred from fluc- 
tuation scale dependence in the left panels by inverting 
Eq. 0} 0| . Autocorrelations have distinct same-side 
(\4>a\ < tt/2) and away-side (|<M| > t/2) components. 
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For peripheral collisions (top-right panel) the same-side 
peak appears to be nearly symmetric on (t)a,<Pa), how- 
ever, cf. the peak widths in Fig. In general, the cor- 
relation structure evolves rapidly with centrality. 

Errors for (p t ) fluctuation measurements in Fig. ^ 
(left panels) are discussed in 0. Statistical errors for 
those variance differences are typically less than 0.001 
(GeV/c) 2 for all scales and centralities. The inversion 
process (effectively a differentiation, which acts as a 
'high-pass' filter) tends to exaggerate small-wavelength 
noise on the autocorrelation. Control of that noise during 
inversion requires a standard procedure called regulariza- 
tion, in which each bin of Ap/ \j p re f is treated as a \ 2 
fitting parameter, incorporating a smoothing term with 
corresponding Lagrange multiplier into the \ 2 expres- 
sion 0, 0| . Autocorrelation errors then have two com- 
ponents: statistical noise which survives smoothing and 
systematic error due to image distortion by smoothing. 
Statistical errors on the autocorrelation are estimated 
by inverting the noise estimate for A<jp t . n . The per-bin 
r.m.s. statistical error which survives smoothing is about 
0.0002 (GeV/c) 2 for all autocorrelations. The smoothing 
distortion, estimated by passing data through the inver- 
sion process twice and c omp aring the resulting two auto- 
correlation versions 0, Il7| . typically peaks at about 5% 
of the maximum autocorrelation value at points of max- 
imum gradient. Correlation amplitudes inferred from 
model fits (see below) were corrected for tracking inef- 
ficiencies and background contamination [lfj| | . An overall 
systematic error of ±14% for corrected amplitudes re- 
flects uncertainty in extrapolation of variance-difference 
measurements to the true number of primary particles in 
the acceptance. 




FIG. 2: Distributions of Ap/^p7^(?7A, (Pa) for 80-90% (up- 
per left), 45-55% (upper right), 20-30% (lower left) and 0- 
5% (lower right) of total cross section. Monopole (constant), 
dipole and quadrupole components have been subtracted. 

In Fig. [21 monopole (constant offset), dipole cos(^a) 



and quadrupole cos(20a) components (sinusoids evident 
in Fig. right panels) have been subtracted from the 
autocorrelations for centralities (80-90, 45-55, 20-30 and 
0-5%) by minimizing residuals of the three components 
on the away side (\4>a\ > 7r /2) and for \t)a\ ~ 2 (mini- 
mizes influence of away-side peak structure). The full-0- 
acceptance fluctuations in Fig. ^ (left panels) are deter- 
mined only by the minijet structure and the monopole 
component; the dipole and quadrupole components inte- 
grate to zero at 5<f) — 2n. Since one interpretation of p t 
correlations is that they reflect velocity correlations of lo- 
cal particle source velocities, the quadrupole component 
from this analysis may constitute the first observation of 
elliptic flow as a true velocity phenomenon. 

The subtracted autocorrelations retain three struc- 
tures localized on both t]a and 4>a- a same-side positive 
peak, a same-side negative peak (apparent as the regions 
of negative correlation immediately adjacent to the posi- 
tive peak on 4>a) and an away-side peak. The near-side 
negative peak cannot be a result of incorrect subtraction 
of the multipole components. The latter have by def- 
inition no structure (are constant) on j]a, whereas the 
negative near-side peak is highly structured (a peak ris- 
ing symmetrically to zero) on 77A ■ The near-side positive 
peak, in the absence of alternative explanations, is inter- 
pretable as a velocity correlation associated with semi- 
hard parton scattering (minijets). Those three pt corre- 
lation structures comprise the main subject of this paper. 

In Fig. |2 we observe that the three peak features vary 
strongly in shape and amplitude with collision central- 
ity. For the more central collisions we observe that the 
same-side positive peak is substantially elongated along 
r/A and significantly narrowed along 4>a- We quantify 
those observations with model fits. The autocorrelations 
in Fig. |2 were fitted with the model function defined in 
Eq. (J2J, a sum of near-side positive peak B\, near-side 
negative peak B2 (signed number) and away-side peak 
B3 , each with the same form, 



F = ex P 



!- 






4>A - S i3 it 


T 4>i 1 













,(2) 



where 8a is a Kronecker delta. This function includes ex- 
ponents t as shape parameters. In contrast to a gaussian 
(r = 2), which best describes near-side peaks for number 
autocorrelations |16| , best-fit exponents for these pt auto- 
correlations were found to be t v i = t, ; 3 = t^i = 1.5 ±0.1, 
with T 2 = 2.5 ± 0.1, T 03 = 1.9 ± 0.1 and t v2 = 1.7 ± 0.1 
(for all centralities). Widths for near-side negative and 
away-side peaks varied (respectively from peripheral to 
central) nearly linearly over the ranges 0.75 < 0^2 < 
1.1 ± 0.1, 0.9 < <r„ 3 < 1 ± 0.1, cr 2 ~ 2.1 ± 0.2 and 
2.4 > ct 03 > 1.5 ±0.1. 

The best-fit amplitudes for all peaks, corrected for 
background contamination and tracking inefficiency (21- 
38% for these 200 GeV data, increasing from minimum 
to maximum centrality) [T(i| and widths for the near-side 
positive peak are plotted in Fig. [3] vs mean participant 
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FIG. 3: Correlation amplitudes (left panel) plotted vs mean 
participant path length v. Solid dots are for positive near-side 
peaks, solid triangles are for negative near-side peaks. Open 
circles are for away-side peaks. Positive-peak widths (right 
panel) are plotted on v for pseudorapidity (solid triangles) 
and azimuth (open circles). Curves guide the eye (see text). 
Dotted, dashed and dash-dot lines represent Hijing results. 



path length v (231 |2J] . The vertical line to the right of 
each panel indicates the estimated limit of v for Au-Au 
collisions, corresponding to b ~ and N part /2 ~ 191. 
The fitted peaks are strongly localized on r]A, and very 
different in shape from the 77A-independent sinusoid com- 
ponents subtracted from autocorrelations in Fig. ^ (right 
panels) to form those in Fig. [3 There is therefore negli- 
gible cross-talk between the two types of fitting function. 
The difference between fitting the data in Fig. ^ (right 
panels) including sinusoids in the model function and the 
data in Fig.|2|with Eq. (J2J is less than the stated errors in 
the fit parameters. Residuals from fitted peak structures 
were typically less than 2% of the near-side peak ampli- 
tude. The peak amplitudes increase with centrality to a 
maximum value and then decrease sharply for the most 
central collisions. The near-side positive-peak width on 
r)A increases monotonically with centrality, while that on 
4>a decreases. The rising part of the B\ data is consistent 
with both a straight line oc (y — 1) and a curve oc v 1G . 

The autocorrelation density ratio Ap(p t : n) / p re f(n) 
measures relative covariances (proportional to Pearson's 
correlation coefficient) of (p t ) fluctuations at pairs of 
points on (j],4>) separated by (t?a>^a)- The autocor- 
relation distribution reveals the average shape of local- 
ized structures which may appear in different places 
on (rj, (j>) in different events, and possibly in only some 
events, but which have some shape stability over an event 
ensemble. We observe substantial (0.01 - 0.1) pt co- 
variances which can be interpreted as local transverse- 
velocity and/or temperature fluctuations correlated at 
pairs of points on (rj, 4>). As with the separation of tem- 
perature and transverse-flow aspects of inclusive single- 
particle p t distributions, temperature and velocity cor- 
relations may also be distinguished as to source mecha- 
nism if mass identification is implemented in pt fluctua- 
tion/correlation analysis. 

The structures in Figs. ^ (right panel) and|21can be 
compared to signature angular correlations for high-pf 
parton fragments: a same-side 2D peak at the origin (jet 



cone) and an away-side ^-invariant ridge on (ap- 
parent as such only in 2D analyses like this one). The 
correlations in Fig. ^ (upper-right panel) have exactly 
that structure, and agree in detail with the Hijing cor- 
relations described below which are known to represent 
\ow-pt jets, providing strong evidence that the dominant 
source of p t correlations for p-p and peripheral Au-Au 
collisions is low-p t parton fragments. For low-pt partons 
the away-side ridge in Fig. (upper-right panel) is not 
distinguishable from one lobe of a sinusoid, and is thus 
removed in the multipole subtraction to obtain Fig. 
The same-side positive peak then represents the conven- 
tional jet-cone structure, albeit in p t rather than angular 
correlations. We can then argue that the same features 
for more central Au-Au collisions continue to derive from 
low-p t partons, but with modifications by a colored dis- 
sipative medium. 

The model-fit results in Fig. [3] illustrate the dramatic 
changes in Fig. [21 structure with collision centrality. Fo- 
cusing on the near-side positive peak, the width on 0a 
falls by 30%, whereas the width on t)a increases by 
60%. Those trends are qualitatively similar to equiv- 
alent measurements of angular correlations |l(| . where 
the large width increase on 7/a was interpreted as due to 
strong coupling of low-pt partons to the longitudinally- 
expanding colored medium. The a v variation is much 
less for p t correlations, suggesting that elongation of par- 
ton fragment distributions on rj involves lower-p t parti- 
cles with increasing rjA- The near-side peak for p t cor- 
relations is significantly non-gaussian , the sharper peak 
represented by exponent r = 1.5 (the correspondent for 
angular correlations is a gaussian with r = 2). 

The amplitude B\ of the near-side positive peak in- 
creases by a factor four or more with centrality, but falls 
off rapidly for the most central collisions, in contrast to 
the subtracted monopole term which increases monoton- 
ically to mid centrality and then remains constant with 
further centrality increase. Since the relative covariance 
could also be interpreted (with strong assumptions) as a 
number of correlated pairs per particle in the system, the 
increase of B\ may indicate that the number of correlated 
pairs from minijets increases relative to the total multi- 
plicity. If the system were composed only of independent 
minijets (with no soft particle production) the autocor- 
relation density ratio (and variance difference) would be 
independent of system size (minijet number). The ob- 
served increase with system size could then result from a 
larger number of low-pt partons, a larger mean fragment 
multiplicity for each parton, or both. The other cor- 
relation structures, the negative same-side and positive 
away-side peaks, are unique to p t correlations and will 
be considered in detail in a followup publication. The 
presence of the negative near-side peak means that the 
variance difference, as an integral fluctuation measure, 
significantly underestimates the relative amount of mini- 
jet structure. 

An equivalent analysis of (pt) fluctuations in Hijing 
collisions ^} indicates that the near-side peak shapes 
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for Hijing (quench on or quench off) are nearly symmet- 
ric on (tja, 0a) j with shape described by single exponent 
r = 1.7 ± 0.1. The combination of same-side 2D peak 
and away-side ID azimuthal ridge observed in that anal- 
ysis supports the interpretation that the basic source of 
those pt correlations is low-p t parton fragments or mini- 
jets, consistent with the basic pQCD jet model in Hijing. 
The centrality dependence of the amplitudes of the Hi- 
jing near-side peak for quench-on (default) and quench- 
off collisions is represented respectively by dotted and 
dashed lines in Fig. (left panel). The lines in the right 
panel correspond to Hijing default (quench on) same-side 
peak widths. The amplitude (width) centrality trends 
for default Hijing are similar: modest variations linear 
with path length v. The differences between quench-on 
and quench-off results from central Au-Au collisions for 
amplitudes and widths, representing pQCD modeling of 
in-medium parton energy loss, are ~ 10%. 

Comparing Hijing to the present analysis in Fig. yjwe 
note that there are at least three instances of qualitative 
disagreement between Hijing and RHIC data. First, for 
very peripheral Au-Au collisions (and therefore nucleon- 
nucleon collisions) the Hijing same-side peak (jet cone) 
is symmetric on angle, whereas the data are quite asym- 
metric. Thus, there is already disagreement with data 
at the level of parton fragmentation in elementary colli- 
sions. Second, the centrality dependencies of the ampli- 
tudes and widths of the same-side peak are qualitatively 
different from data: in some cases even the sign of the 
variation is wrong. Third, the amplitude of the same- 
side Hijing peak is qualitatively less than that for data 
in mid-central collisions. The last is especially surprising 
when comparing quench-off Hijing (dashed line in the left 
panel) to data. The quench-off Hijing option in principle 
models jet production from a linear superposition of N-N 
collisions combined with a Glauber model of a transpar- 
ent nucleus. That model should provide an upper limit 
for jet structure in A-A collisions. Yet the same-side 
peak amplitude for quench-off Hijing is 2-3 times less 



than that for mid-central RHIC Au-Au collisions which 
are observed to be highly opaque to minijets in the central 
region. Finally, we observe no evidence in Hijing data for 
the same-side negative peak which is a prominent new 
feature of RHIC data. The perturbative treatment of 
parton energy loss in Hijing appears to disagree strongly 
with the observed process for the minimum-bias partons 
which dominate pt correlations in Au-Au collisions. 

In conclusion, we have for the first time measured the 
scale dependence of (pt) fluctuations on {5rj, 5<t>) in heavy 
ion collisions. We have inverted those distributions to ob- 
tain autocorrelation distributions on corresponding dif- 
ference variables (^Aj^a) which reveal the correlation 
structure of the local properties of a two-particle pt dis- 
tribution, specifically a combination of local transverse 
velocity and temperature. Inferred autocorrelations re- 
veal complex p t correlation structure in Au-Au collisions 
at RHIC, including peaked structures attributed to mini- 
jets which vary strongly with collision centrality. We 
observe that pt autocorrelations provide unique access 
to minijet structure down to very low pt and probe the 
detailed interplay between \ow-p t partons and the dissi- 
pative colored medium. Further studies with identified 
particles may separately characterize the local velocity 
and temperature structures of heavy ion collisions. 
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